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ABSTRACT 

We present stellar metallicities derived from Ca II triplet spectroscopy in 
over 350 red giant branch stars in 13 fields distributed in different positions in 
the SMC, ranging from ~1° to ~4° from its center. In the innermost fields the 
average metallicity is [Fe/H] ~ —1. This value decreases when we move away 
towards outermost regions. This is the first detection of a metallicity gradient in 
this galaxy. We show that the metallicity gradient is related to an age gradient, 
in the sense that more metal-rich stars, which are also younger, are concentrated 
in the central regions of the galaxy. 
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Introduction 



The chemical enrichment history of a galaxy is related to the origin and distribution of 
nuclear species in its stars and gas. The chemical elements are mainly produced by stars 
which drive the enrichment of the interstellar medium by ejecting material containing the 
product of the stellar nucleosynthesis from which the new generations of stars are formed. In 
addition, gas flows also play an important role in chemical enrichment, diluting the products 
of the stellar nucleosynthesis with unenriched material from outside the galaxy, and mixing 
metals from one part of the system to another (e.g. bringing metal-rich gas into metal- 
poor regions). Thus, the study of chemical evolution of galaxies involves understanding the 
spatial distribution and temporal evolution of various elements by taking into account the 
processes of star formation, the distribution of stars according to their masses and chemical 
compositions, and the final yields of various elements and detectable remnants of parent 
stars. Until recently, only the chemical enrichment history of the solar vicinity could be 
studied in detail. However, the modern multiobject spectrographers attached to the 8-10 
m class telescopes allow us to study the chemical enrichment history of the nearest Local 
Group galaxies. 

Because of their proximity, and the fact that they present a wide range of ages and 
metallicities, the Magellanic Clouds are attractiv e objects to study chemical enrichment 
histories. In a previous paper (jCarrera et al.ll2008l . hereafter paper III), we investigated the 
chemical enrichment history of the LMC. In this paper, we will focus in the study of the SMC. 
There are considerable less studies of the SMC as compared with the LMC, probably due to: 
(i) its irregular appearance, with complex kinematics; (ii) its distance, located further away 
than the LMC; and (iii) its depth in the line of sight, which remains a subject of controversy. 

Most of the inform ation about the stellar population s of the SMC has been ob tained from 



its cluster system (e.g. iDa Costa &: Hatzidimitrioulll998l : iPiatti et al 



ter age distribution does not show the age- gap observed in the LMC ( jPa Costa fc Hatzidimitriou 



1998 



2001 



2005). The clus- 



Rich et al 



Mighell. Saraiedini. fc Frenchlll998l ). From a sample of seven clusters older than 1 Gyr, 



(120001 ) suggest that star formation was stronger in two main episodes, one ~8±2 
Gyr ago and another ~2±0.5 Gyr ago. However, there does not seem to be any age interval 
completely lacking of objects (IRafelski fc Zaritskyll2005l ). Ther e is only one old cl uster, NGC 
121, that is younger than most of the LMC globu lar clusters (IPiatti et al.ll2005l ). However, 



the SMC shows a significant old field population (INoel et al.l 120071 . hereafter Paper I) 



To our knowledge, there are only thr ee studies in which the star formation history (SFH) 



of the SMC field population was derived ( Dolphin et al. 
a deep color-magnitude diagram (CMD) iDolphin et al. 



20011: 



Harris fc Zaritskyll2004l ). From 



(120011 ) found than the star formation 



in a small field in the periphery of this galaxy was relatively constant until about 2 Gyr ago. 
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with no star formation occurring since then. However, this could be bias ed by the fact that 
their field was specifically chosen because it did not have young stars. [Harris fc Zaritsky 
( 120041 ) studied a 4°x4?5 field in the central region of the galaxy. From shallower CMDs they 
found that the star formation in the SMC has had two main episodes: one which formed 
the oldest populations and lasted until 8.5 Gyr ago, and a recent one that started around 3 
Gyr ago. Noel et al (2008, in preparation) have also obtained accurate SFHs for the fields 
presented in this paper, using CMDs reaching the oldest main sequence turnoffs with good 
photometric accuracy. They find that two main episodes of star formation in all fields, one 
at old ages (~ 10 Gyr ago) and another one at intermediate ages (~ 5 Gyr ago), in addition 
to young star formation in the wing fields. 

Detailed det erminations of chemical abundances exist only for the youn gest population 



of the SMC (i.e. Iffill. Barbuv. fc Spitd 119971 : IVennlll999l : iHunter et al.ll2007f ). and its chem- 



ical enrichment history has been mainly determined from studies from its cluster systern . 
The cluster age-met allicity relationship (AMR) has been obtained by iPiatti et al.l (120051 ) 
and iMighell. Sarajedini. fc FrenchI (119981 ). mainly from photometric indicators. An initial 
chemical enrichment has been found, followed by a period of relatively slow increase in the 
metal abundance. Clusters more metal-rich than [Fe/H]>-1 are younger than 5 Gyr. Since 
then, the metallicity has again increased until now. On average, t he SM C is more metal-poor 
than the LMC. The very recent work by lldiart. Maciel fc Costal (120071 ). which have obtained 
chemical abundances in a sample of SMC planetary nebulae, has found a similar result, with 
the exception that the chemical enrichment episode at a very early epoch is not observed. 

In the present work we focus on obtaining stellar metallicities of individual red giant 
branch (RGB) stars in the field population of the SMC from Ca II triplet (hereafter CaT) 
spectroscopy. These stars have been selected in 13 fields distributed at different positions in 
the SMC ranging from ~1° to ~4° from its center. Deep photometry of these fields has been 
presented in Paper I . The procedure followed to select the targets is explained in Section [2J 
The data reduction is discussed in Section [31 The radial velocities of the stars in our sample 
are obtained in Section HI In Section [5] we discuss the calculation of the CaT equivalent 
widths and the determination of metallicities. Section [6] presents the method used to derive 
ages for each star by combining the information on their metallicity and position on the 
CMD. The analysis of the data is presented in Section [71 The metallicity distribution of 
each field and the possible presence of a metallicity gradient is discussed in Section 17.11 The 
derived AMRs for each field are presented in Section 17.21 The main results of this paper are 
discussed in Section [SI 
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2. Target Selection 

In the framework of a large program to obtain proper motions, deep CMDs and stel- 
lar metallicities in the SMC, we secured spectroscopy of stars in 13 fields spread about 
the galaxy body. The photometry of these fields, presented in Paper I and listed in Ta- 
ble [H were obtained with a Tektronic 2048x2048 CCD detector attached to the LCO 100" 
telescope, which covers a field size of 8f85x 8^85. Following the notation described by 



Tinney. Da Costa &: Zinneckerl (jl997l ). fields denoted by "qj" (followed by right ascension) 
are centered on quasars and were observed photometrically with the main objective of deter- 
mining the absolute proper motion of the SMC (Costa et al., in preparation). Fields labeled 
"smc" were selected specifically to study their stellar populations by sampling a range of 
galactocentric radius at similar azimuth. The BR photometry of these fields is described 
in detail in Paper I, where the distribution of stellar populations of the SMC is discussed 
on the basis of CMDs. These observations have been complemented with observations in 
the I band in order to allow using t h e red uced equivalent width- met allicity (W'/-[Fe/H]) 



relationship derived in ICarrera et al.l (120071 . hereafter Paper 11) to obtain metallicities for 
individual RGB stars, /-band observations of "qj" fields were obtained with the same in- 
strument and telescope as the BR photometry. The / images of "smc" fields were obtained 
with F0RS2 at the VLT in image mode and were also used for spectroscopic mask configura- 
tion of MXU@F0RS2. Field qjOlll was observed with both telescopes in order to compare 
the photometric calibrations. The magnitudes obtained with each calibration differ by about 
~0.1 magnitudes. This relatively large difference is mainly explained by the poor photomet- 
ric calibration of the F0RS2 images, which were not taken for photometric purposes. In any 
case, this error means a metallicity uncertainty of only ~0.02 dex, a value smaller than the 
metallicity uncertainty itself (~0.1 dex). 

In each field we selected stars in two windows of the CMD, which are plotted in Figure 
[U In each region the stars were ordered from the brightest to the faintest ones, regardless 
of their color. The resulting star list was used as input for the mask configuration task of 
the instrument. Stars in the box below the RGB tip were given higher priority, and only 
objects above the tip were observed when it was impossible to put the slit on a star of the 
lower region. 



3. Observations and Data Reduction 

The spectroscopic observations were carried on in service mode with the VLT ANTU 
telescope, at Paranal Observatory (Chile), through program 074.B-0446. We used F0RS2 
in MXU multiobject mode with grism 1028z+29 and filter OG590+32 in order to eliminate 
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residual orders. Slits with a length of 8" and a width of 1" were selected. With this setup, we 
were able to observe about 40 objects in each field (the final number depends on the spatial 
distribution of the selected stars). The slit length was selected with the purpose of shifting 
the stars along it in order to acquire two exposures of each field without superposition of 
the stellar spectra in each. This particular procedure allowed to extract the spectra in the 
same way as Paper II. 

After bias subtraction, each image was fiat-field corrected. Then, as each star is in a 
different position in the two images of the same field, we subtracted one from the other, 
obtaining a positive and a negative spectrum of the same star. With this procedure the 
sky is subtracted in the same pixel in which the star has been observed, thus, minimizing 
the effects of pixel-to-pixel sensitivity variations. Sky residuals due to temporal variation of 
the sky brightness were eliminated in the following step, in which the spectrum is extracted 
in the usual way and the remaining sky background is subtracted using the information on 
both sides of the stellar spectrum. In the next step, the spectra were wavelength calibrated 
and added to obtain the final spectrum. Finally, each spectrum was normalized by fitting a 
polynomial, excluding the strongest lines (such as the CaT lines). There is an uncertainty 
in the wavelength calibration because the arcs used for this purpose were not taken at 
the same time and with the same telescope pointi ng as the obje c t. Th e effects of this on 



the wavelength calibration have been discussed by iGallart et all (l200ll ). Since we are not 
interested in obtaining precise radial velocities, this problem will not significantly affect our 
results. 

In total, we observed 386 stars in the 13 SMC fields. Their magnitudes and CaT 
equivalent widths are listed in Table [2] 



4. Radial Velocities 

The radial velocity of each star was calculated in order to reject SMC non-members. 
We used the IRAF fxcor task, which performs the cross-correlation betw een the target and 



template spectra of stars of known radial velocity (ITonry fc Davislll979l ). As templates we 
selected nine stars in the clusters NGC 104, NGC 2682, NGC 288 and NGC 7078 which 
were observed within the same program as the SMC fields and presented in Paper II. The 
velocities were corrected to the heliocentric reference frame within fxcor. The final radial 
velocity for each SMC star is the average of the velocities obtained from each template, 
weighted by the width of the corresponding correlation peaks. The resulting velocities can 
be affected by the fact that the stars might have not been positioned exactly in the center 
of the slit. The importance of this is described in detail in Paper II. However, in the case 
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of the SMC, the stars were observed at a relatively large air mass (>1.6) and with a seeing 
near to or larger than 1". Since the slit width was 1", the effect of the incorrect centering of 
the star on the slit was of little importance. In fact, when we tried to characterize this effect 
in the way described in Paper II, we found that its contribution to the resulting velocity is 
smaller than the uncertainty due to the wavelength calibration. For this reason, we did not 
take this effect into account in the final radial velocity. We cons idered as SMC member s 



those stars with radial velocities in the range 50<Vr <250 km s ^ (IHarris fc Zaritskyll2006l ). 



5. CaT Equivalent Widths and Metallicity determination 

The metallicity of the RGB stars is obtained following the procedure described in Paper 
II. The equivalent width is the area of the line normalized to the local continuum within 
a line bandpass. The continuum is calculated from a linear fit to the mean value of the 
corresponding bandpasses. The line and continuum bandpasses used in this work are listed 
in Table [31 The line flux is calculated from the fit of its profile using a Gaussian plus a 
Lorentzian function. As discussed in Paper II, this function provides a better fit to the 
core and the wings of the strongest lines than other functions previously used. The CaT 
index, denoted as ECa, is defined as the sum of the equivalent widths of the three CaT 
lines. The HCa and their uncertainties for each star observed are given in Table [21 together 
with their magnitudes and radial velocities. Two calibrations of the CaT as metallicity 
indicator were obtained in Paper II based on / and V magnitudes. In this case, only / 
magnitudes are available for the SMC stars. The reduced equivalent width, W'/, for each 
star has been calculated using the slope obtained in Paper II for the calibration clusters in 
the Mj-SCa plane (/3/=-0.611 Amag~^). To obtain the abso lute magnitudes we assumed 



a distance modulus of (m-M)o=18.9 (see Ivan den Berghlll999l ) and reddenings listed in the 



last column of Table [H (see Paper I for details). Also, in Paper II three different metallicity 
scales were used as r eferen ce. In this case, we only used the relationships obtained on the 



Carretta fc GrattonI (Il997l . hereafter CG97) metallicity scale, because it is the only one that 
uses homogeneous high-resolution metallicities of open and globular clusters, and because 
the metallicities of the LMC stars in Paper III were also obtained in this way. 

In brief, the metallicity for each star is given by: 

[Fe/H]cG97 = -2.95 + 0.38SCa + 0.23M/ (1) 



In Figures [21 [31 and [H the position of SMC stars (radial velocity members) in the M/- 
SCa plane for our eastern, western, and southern fields, are respectively shown. Solid lines 
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indicate the magnitude range of the observed cluster stars used in Paper II to obtain the 
above relationship. Dashed lines show the region where the relationship is extrapolated. 

The metallicity distribution of each field is shown in Figures El M and [7] and will be 
discussed in Section 17. 1[ 



6. Determination of stellar ages 

The position of the RGB on the CMD suffers from age-metallicity degeneracy. However, 
when the metallicity is obtained in an alternative way, as in this case from spectroscopy, 
this age-metallicity degeneracy can be broken, and ages can be derived from the position 
of the stars in the CMD. In paper III, a polynomial relationship was computed to derive 
stellar ages from their metallicitie s and positions in the CM D. For that purpose, synthetic 



CMDs computed with lAC-star (lAparicio fc GallartI 120041) with the overshooting BaSTI 



(jPietrinferni et al.l |200J) and Padova (jCirardi et al.l |2002| ) stellar stellar evolution models 



as input were used. As it is explained in Paper III, differences between both models in 
the resulting ages are negligible for our purpose, since we are not interested in an accurate 
determination of ages. For simplicity, as in Paper III, we used only the relationship obtained 
from the BaSTI stellar models. The aforementioned relation was obtained for {V — I) and 
My. Since V magnitudes are unavailable for our sample of SMC stars, we computed a new 
relationship for B, R, and I magnitudes. In this case, we selected {B — I) instead of {B — R) 
because the former is much more sensitive to small changes in the stellar metallicity and age. 

We followed the same procedure as in Paper III. First, we used the same synthetic CMD 
used in Paper IIl[^, which was computed with a constant star formation rate (SFR) between 
and 13 Gyr and with a chemical law such that any star can have any metallicity between 
—2.3 to +0.5 dex. In this CMD we only selected stars in the same region below the tip of the 
RGB in which the observed ones were mainly chosen. We did not consider the brightest AGB 
stars due to the uncertainty of their parameters as predicted by stellar evolution models. 
Following the same statistical procedure described in Paper III, we checked which polynomial 
combinations of magnitude Mj, color {B — J) and metallicity [Fe/H] best represent the age 



^During the rcfcrrcing process of this paper, the authors of the BaSTI models discovered a problem in 
the calculation of models in the mass range 1.1-2.5 Mq (~1-4 Gyr), (see littp:/ /albione.oa-teramo.inaf. it/). 
Even though the differences between the new and old models were unlikely to affect in any substantial way 
our results, we recalculated equation [2] using a synthetic CMD computed with the updated model set. All 
the results shown are derive from this new relationship. We have verified that the differences are indeed 
minor, and that, therefore, the results on the LMC on paper III, obtained using the faulty models, can be 
trusted nevertheless. 
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of the stars in this synthetic CMD. In order to minimize the a and to improve the correlation 
coefficient, different hnear, quadratic and cubic terms of each observed magnitude have been 
added. We have checked whether M/ or Mb magnitudes improve the relationship. Similar 
results are obtained for both magnitudes, so we choose the ffist because metallicities were 
calculated from it. The final polynomial form adopted is: 

log{age) = a + b{B-I) + cMj + d[Fe/H] + f{B - if + g[Fe/HY (2) 
The best fit coefficients are listed in Table HI 

In order to estimate the age uncertainty when this relationship is used to compute stellar 
ages, we performed a Monte Carlo test as in paper III. The goal is to check how the obtained 
ages change when the input parameters are modified. The test consists in computing, for 
each synthetic star, several age values for stochastically varying [Fe/H], [B — /) and M/ 
according to a gaussian probability distribution of the corresponding a {criFe/H] ~0.15 dex; 
^{B-i) ~0.001 and ~ 0.001). The a value of the obtained ages provide an estimation of 
the age error when Equation [2] is used. The obtained values for the considered age intervals 
are shown in Figure [HI The age uncertainty increases for older ages. 

It is also necessary to check how well equation [2] reproduces the age of a real stellar 
system. Following the same steps as in Paper III, we choose the cluster stars used in Paper 
II for the calibration of the CaT as metallicity indicator. For those stars we knew their 
[B — I) color and Mj magnitudes, so we can compute their metallicity in the same way as 
was done for the SMC stars (see Section [5]). We then used these observational magnitudes as 
input for equation [2|, obtaining an age for each cluster star. The cluster age was computed 
as the mean of the ages of its member stars. In Figure [HI the age computed for each cluster 
has been plotted versus its reference value. As in Paper III, ages younger than 10 Gyr are 
well recovered. However, the relationship saturates for ages larger than 10 Gyr. 

7. Analysis 

7.1. Metallicity distribution 

Metallicity distributions are shown in Figures O [6] and [7] for eastern, western and south- 
ern SMC regions respectively. We have ffited a Gaussian to obtain the mean metallicity and 
metallicity dispersion of each of them. These values are listed in columns 3 and 4 of Table 
[Si The fields are ordered by its distance to the center, which is shown in column 2. Fields in 
different regions are indicated by different font types: eastern fields in normal, western fields 
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in boldface and southern fields in italics. Mean metallicities are very close to [Fe/H]~-1 in all 
fields within r<2?5 from the SMC center. A similar value is observed for the southern fields 
up to r<3° (qj0047 and smc0049). Note that the SMC isopleths of intermediate-age and 



old stars are elongated in the NE-SW direction (PA=45°; ICioni. Habing fc Israeli |2000| ) and 
that a radius of 3° in the southern direction correspond to approximately the same isopleth 
at radius 2?5 in the eastern and western directions. For the outermost fields, qj0033 in the 
West, and qj0102 and qj0053 in the South, the mean value is clearly more metal-poor than 
in the others. 

The fact that the mean metallicity decreases when we move away from the center implies 
that there is a metallicity gradient in the SMC. This gradient is more clear when we compute 
the percentage of stars in different metallicity bins, values also listed in Table O (columns 5 
and 6). For the western and southern fields, the percentage of stars more metal-poor than 
[Fe/H] = —1 increases when we move away from the center. This is not observed in the 
eastern fields because they are almost at the same distance. This is the first time that a 
metallicity gradient has been reported in SMC stellar populations. The detection of this 
gradient has been possible because we have covered a large radius range, up to 4° from the 
SMC center. 



7.2. Age— metallicity relationships 



From a purely phenomenological point of view, there are two main ways to account 
for the mean metallicity gradient found in the previous section. One possibility is that 
chemical enrichment has proceeded more slowly towards the SMC periphery, in such a way 
that coeval stars would be more metal-poor when we move away from the center. This seems 
to be the case of spiral galaxies, like the Milky Way, where the observed abundance gradients 
may be explained by radial variations of the relation between the SFR and the am ount of 
infalling gas (e.g. iPrantzos Sz Boissierl l2000l : IChiappini. Matteucci Sz Romano! 1200 ll ). The 



situation would be complicated in dwarf galaxies by the probable existenc e of galactic winds 
originated in supernova explosions (e.g. iRomano. Tosi &: Matteucci! l2006l ) which are able to 
remove large amount of metals from the interstellar medium. An alternative scenario is that 
the stellar AMR (i.e. the law of chemical enrichment as a function of time) has been the 
same everywhere in the SMC. As a result, the average metallicity of coeval stars would be 
the same in all fields and the metallicity gradient would be related to an age gradient. A 
mixture of both scenarios is also possible. 



To investigate the nature of the gradient, we have therefore calculated the AMR for 
each field. They are plotted in Figures !TU l !TT | and !T^ for fields situated to the East, West 
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and South respectively. Note that the uncertainty in age is much larger than in metallicity. 
However, as we are interested in the global behavior and not in obtaining individual stellar 
ages, the age determinations are still valid. Since the procedure to obtain the age saturates 
for values older than 10 Gyr, we can only be confident that the oldest stars have an age > 
10 Gyr. For these we as sume an age of 12.9 Gy r, which is the age of the oldest cluster in the 



Milky Way (NGC 6426. ISalaris &: Weissll2002l ). Regarding the youngest stars, in the region 



of the CMD where we selected the observed stars, and according to the stellar evolution 
models, we do not expect to find stars younger than ~0.8 Gyr. However, equation [2] can 
formally compute ages younger than this value. As the age determination uncertainty for 
these young stars is ~1 Gyr, in order to avoid this contradiction we assign them an age 
of 0.8 Gyr. Inset panels show the age distribution of the observed stars, with and without 
taking into account the age uncertainty [solid line and histogram, respectively). To obtain 
the first one, we assumed that the age of each star is represented by a Gaussian probability 
distribution on the age axis, with a mean value equal to the age calculated for this given 
star, and a equal to the age uncertainty. The area of each of these distributions is unity. In 
the case of stars near the edges, the wings of the distribution may extend further than the 
limits. We have proceeded in the same way as described in Paper HI, cutting off the wings 
outside the assumed limits (0.8 and 12.9 Gyr) and rescaling the rest of the distribution so 
that the area remains unity. 

All the AMR plotted in Figures [TOl [TT] and [T2] show a rapid chemical enrichment at a very 
early epoch. Even though in some fields we have not observed enough old stars to sample this 
part of the AMR, note that 12 Gyr ago all fields have reached [Fe/H]~ —1.4 — — 1.0. This 
initial chemical enrichment was followed by a period of very slow metallicity evolution until 
around 3 Gyr ago. Then, the galaxy started another period of chemical enrichment, which is 
observed in the innermost fields, which are, however,the only ones where we observed enough 
young stars to sample this part of the AMR. For a given age, the mean metallicity of the 
stars in fields qjOlll, qj0112 and smc0057 seems to be slightly more metal-poor than those 
of other fields. The uncertainties in the age determination could account for the observed 
differences. In all cases, however, the mean metallicity is similar to that of the other fields 
at similar galactocentric radii. The field AMRs obtained in this work are similar to those for 
clusters (the reader should take into account that differences in the metallicity scales exist 
among different works), although there is only one cluster older than 10 Gyr (see Figure [T3|l . 
and for planetary nebulae, with the exception that in these object s it is not observed the 



chemical enrichment episode at a very early epoch (see Figure 6 of lldiart. Maciel fc Costa 



20071 ). 



It does not seem that there was a period in which the galaxy has not formed stars, in 
agreement with the result found in Paper I. (Figures \TU \ [TT] and [T2 | inset panels). For eastern 
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fields, located in the wing, most of the observed stars have ages younger than 8 Gyr, but 
there is also a significant number of objects older than 10 Gyr. At a given galactocentric 
distance, eastern fields show a large number of young stars (<3 Gyr) in comparison to the 
western ones, as discussed in Paper I. For the western and southern fields, the fraction of 
intermediate-age stars, which are also more metal-rich, decreases as we move away from 
the center, although the average metallicity in each age bin is similar. This indicates the 
presence of an age gradient in the galaxy, which may be the origin of the metallicity one. It 
is noticeable that for the most external fields, qj0033 and qj0053, we find a predominantly 
old and metal-poor stellar population. 

We can check statistically the hypothesis that the AMR is independent of the position 
in the SMC. To do so, we have combined the measurements on the 13 fields to obtain a 
global AMR and compared it with that of each field. To do so, we have divided the age 
range into six intervals (age<1.5 Gyr, 1.5-3.5 Gyr, 3.5-5.5 Gyr, 5.5-8.5 Gyr, 8.5-11 Gyr, 
>11 Gyr). We have computed the dispersion and mean metallicity in each age bin, and 
listed the results in Table [61 With these data, we have performed a test as follows: 



where af is the squared sum of the uncertainties in the age bin i of each field and the 
global AMR. The result is shown in Column 8 of Table [61 All fields, except qj0047 and 
smc0033, have values of xl smaller than 1. The discrepancy in fields qj0047 and smc0033 
may be due to small number statistics: there are only two stars younger than 11 Gyr in field 
qj0047 and three in smc0033. Values of xl <1 niean that the observed AMR for each field 
is the same as the global one with a confidence of 90% (95% in most cases). From these 
results, we may conclude that the hypothesis that the AMR is independent of position is 
correct within the uncertainty. 



Using CaT spectroscopy, we have derived stellar metallicities for a large sample of RGB 
field stars in 13 regions of the SMC situated at different galactocentric distances and positions 
angles. We have found a radial metallicity gradient, which is most evident for those fields 
situated toward the South, where we covered a large galactocentric radius. For a given 
galactocentric distance, the mean metallicities for fields situated at different position angles 
are very similar. The inner fields have a mean metallicity of [Fe/H] ~ —1, which is similar 




(3) 



8. Summary and Discussion 
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to that of the cluster metalhcity distribution. 

We have obtained the AMR of each field from the combination of metallicities, derived 
from CaT spectroscopy, and the position of stars in the CM P. All fields have similar AMRs, 



which are also similar to the clusters one (jPiatti et al.ll2005l ). All show a rapid initial increase 



of metallicity, followed by a very slow chemical enrichment period. A new relatively fast 
chemical enrichment episode is observed in the last few Gyrs in the fields within ~2° of 
the center with enough young stars to sample it. From the information on the AMRs, we 
conclude that coeval stars have the same metallicity everywhere in the SMC. The observed 
metallicity gradient is therefore related to an age gradient, because the youngest stars, which 
are also the most metal-rich, are concentrated in the central regions of the galaxy. 

In a forthcoming paper we will try to reproduce the observed AMR with chemical 
evolution models using accurate SFRs, as a function of time, which are being derived by our 
group in each field (Noel et al. 2008, in preparation). 
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Fig. 1. — Color-magnitude diagram of the field qj0035 showing the RGB windows used to 
select the candidates to be observed spectroscopically. 
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Fig. 2. — Position in the M/-SCa plane of observed SMC stars in the eastern fields. Only 
stars with confirmed membership from their radial velocity are represented. The typical 
SCa error is shown in the bottom right corner of each panel. Isometallicity lines, obtained 
from the relationship based on Mj presented in Paper II, have been plotted for reference. 
The solid part of the line is the magnitude interval covered by the cluster stars used for 
the calibration (see Paper II). The dashed part is the region in which the calibration is 
extrapolated. Distances from the SMC optical center are given in the bottom right corner. 
The innermost field is on the left and the outermost one is on the right. 
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Fig. 3. — The same as Figure [21 for the Western fields. They are ordered from the innermost 
field (top left), to the outermost one {bottom right). 
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Fig. 4. — The same as Figure O for the Southern fields. They are ordered from the innermost 
field (top left), to the outermost one {bottom right). 
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Fig. 5. — Metallicity distributions for the eastern fields in our sample. A Gaussian has been 
fitted to each distribution in order to obtain its mean and dispersion. The values obtained 
are shown in each panel. 
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Fig. 6. — The same as Figure [5] for the western fields. 
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Fig. 7. — The same as Figure [5] for the southern fields. 
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Fig. 8. — Age uncertainty as a function of age as calculated through a Monte Carlo test and 
using equation [2l See text for details. 
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Fig. 9. — Ages derived from equation [2] for the cluster sample presented in Paper II, plotted 
against the reference values. The solid line corresponds to the one-to-one relation. 
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Fig. 10. — Age-metallicity relationships for the eastern SMC fields in our sample. Inset 
panels show the age distribution computed taking {solid line) and not taking [histogram) 
into account the age determination uncertainties. Top panels show the age error in each age 
interval (see Figure [H]). Right panel shows the metallicity uncertainty in each metallicity bin. 
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Fig. 11. — The same as Figure [TOl for the western fields. 
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Fig. 12. — The same as Figure fTOl for the southern fields. 
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Fig. 13.— AMR for SMC clust ers. 
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( 120041 . HZ04) in the central region of the galaxy. The mean metallicity in six age bins of our 
global SMC AMR has also been plotted [open circles). Note that the metallicity scales of 
each work may not be exactly the same. 
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Table 1. SMC observed fields 



Field 


02000 


52000 


r(') 


PA (°) 


Zone 


E(B-V) 


smc0057 


00:57 


-73:53 


65.7 


164.4 


South 


0.09 


qj0037 


00:37 


-72:18 


78.5 


294.0 


West 


0.07 


qjoose 


00:36 


-72:25 


79.8 


288.0 


West 


0.07 


gjOlll 


01:11 


-72:49 


S0.9 


S9.5 


East 


0.09 


gj0112 


(71:12 


-72:36 


S7.4 


Sl.O 


East 


0.09 


qj0035 


00:35 


-72:01 


95.5 


300.6 


West 


0.05 


gjOlie 


01:16 


-72:59 


102.5 


95.2 


East 


0.08 


smcOlOO 


01:00 


-74:57 


130.4 


167.5 


South 


0.05 


qj0047 


00:47 


-75:30 


161.7 


187.7 


South 


0.05 


qjOOSS 


00:33 


-70:28 


172.9 


325.0 


West 


0.03 


smc0049 


00:49 


-75:44 


174.8 


184.6 


South 


0.06 


qj0102 


01:02 


-74:46 


179.5 


169.4 


South 


0.05 


smc0053 


00:53 


-76:46 


236.3 


179.4 


South 


0.06 
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Table 2. Red giants stars observed 



«2000 


<52000 


E Ca 




B 


I 


Vr{kms ^) 


Comments 




00:33:53.8 


-70:24:31.3 


2.75 


0.23 


20.26 


18.45 


-307.1 


68.5 


No member 


00:34:00.0 


-70:28:22.9 


2.56 


0.42 


20.31 


18.57 


79.0 


4.2 




00:34:03.1 


-70:27:07.7 


1.89 


0.71 


20.82 


18.97 


78.9 


3.3 




00:34:25.9 


-70:27:47.1 


5.64 


0.19 


19.66 


17.61 


133.5 


3.7 




00:34:06.4 


-70:27:39.3 


1.95 


0.12 


18.67 


17.81 


162.2 


2.9 




00:33:37.1 


-70:27:31.8 


5.10 


0.08 


18.82 


15.82 


69.9 


3.1 




00:33:41.8 


-70:27:14.3 


5.17 


0.20 


19.72 


17.46 


31.7 


2.9 


No member 



Note. — Table [2] is published in its entirety in the electronic edition of Astronomical Journal. A portion 
is shown here for guidance regarding its form and content 



Table 3. Line and continuum bandpasses 



Line Bandpasses (A) 


Continuum bandpasses (A) 


8484-8513 


8474-8484 


8522-8562 


8563-8577 


8642-8682 


8619-8642 




8799-8725 




8776-8792 



Table 4. Coefficient s of the fit of equation [2] for BaSTI stellar evolution models 



(jPietrinferni et al.ll2004l ). The fit standard deviation of the fit (a) is shown in the last 

column. 



a b c d c fa 

4.06±0.03 3.58±0.02 0.663±0.003 -1.314±0.006 -0.402±0.002 -0.077±0.005 0.40 
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Table 5. Metallicity distribution in each field. 



Field 


r(') 


{[Fe/H]) 




[Fc/H]<-1 


[Fc/H]>-1.0 


smc0057 


1.1 


-1.01 


0.33 


53 


47 


qj0037 


1.3 


-0.95 


0.17 


35 


65 


qjoose 


1.3 


-0.98 


0.25 


49 


51 


gjOlll 


1.3 


-1.08 


0.21 


64 


36 


gj0112 


1.4 


-1.16 


0.32 


69 


31 


qj0035 


1.6 


-1.09 


0.24 


66 


34 


gjoue 


1.7 


-0.96 


0.26 


43 


57 


smcOlOO 


2.2 


-1.07 


0.28 


61 


39 


qj0047 


2.7 


-1.15 


0.27 


68 


32 


qj0033 


2.9 


-1.58 


0.57 


85 


15 


smc0049 


2.9 


-1.00 


0.28 


53 


47 


qj0102 


3.0 


-1.29 


0.42 


74 


26 


smc0053 


3.9 


-1.64 


0.50 


92 


8 



Table 6. Average metallicity in six age bins. Also listed are the values obtained from the 
combination of the 13 fields. The last column shows the value xl obtained from the 
comparison of each field with the global one. 



Field 


{[Fe/H]<i s) 


(lFe/H]i,5 


~3.5> 




-5.5> 




-S.5> 


<[f'<!/-Hls.5-ll> 


([Fe/H]>ii) 




smc()()57 


-0.68±0.11 


-0 


98±0 


20 


-1 


09±0 


14 


-1 


34±0 


04 


-1 


24±0.05 


-1.52±0 


24 


0.39 


qj0037 


-0.15±0.22 


-0 


83±0 


10 


-0 


89±0 


17 


-1 


06±0 


11 


-0 


89±0.02 


-1.53±0 


55 


0.44 


qj0036 




-0 


76±0 


07 


-0 


82±0 


17 


-1 


01±0 


12 


-0 


96±0.16 


-1.32±0 


23 


0.17 


gjOlll 


-0.55±0.29 


-0 


96±0 


19 


-1 


07±0 


07 


-1 


17±0 


13 


-1 


56±0.04 


-1.98±0 


22 


0.62 


gj0112 


-0.85±0.34 


-1 


03±0 


28 


-1 


23±0 


06 


-1 


38±0 


20 






-1.77±0 


27 


0.76 


qj0035 




-0 


66±0 


24 


-0 


73±0 


26 


-0 


98±0 


19 


-0 


94±0.14 


-1.27±0 


26 


0.47 


gjOlie 


-0.71±0.68 


-0 


56±0 


12 


-0 


90±0 


15 


-1 


06±0 


20 


-1 


03±0.07 


-1.22±0 


26 


0.30 


smcOlOO 




-0 


80±0 


20 


-0 


91±0 


16 


-0 


89±0 


02 


-1 


16±0.06 


-1.30±0 


26 


0.24 


qj004T 




-0 


63±0 


01 


-0 


67±0 


01 












-1.20±0 


38 


1.07 


qj0033 
















-0 


74±0 


01 


-1 


05±0.18 


-1.75±0 


40 


1.07 


smc0049 




-0 


74±0 


14 


-0 


87±0 


06 


-0 


92±0 


02 


-1 


13±0.04 


-1.28±0 


17 


0.24 


qj0102 




-1 


00±0 


28 


-0 


99±0 


18 


-1 


06±0 


38 


-1 


31±0.15 


-1.49±0 


31 


0.17 


smc0053 






















-1 


27±0.05 


-1.67±0 


43 


0.36 


Global 


-0.68±0.37 


-0 


86±0 


22 


-0 


94±0 


20 


-1 


10±0 


21 


-1 


11±0.21 


-1.42±0 


37 





